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, The lepton asymmetry created in the out-of-equilibrium decay of a heavy 

Major ana neutrino can generate the cosmological baryon asymmetry Yb when 

processed through fast anomalous electroweak reactions. In this work I examine 

Q_i| this process under the following assumptions: (1) maximal Vn-v T mixing (2) 

hierarchical mass spectrum m% ~ 5 x 10~ 2 eV S> mi (3) small-angle MSW 

solution to the solar neutrino deficit. Working in a basis where the charged lepton 
• , 

^ ■ and heavy neutrino mass matrices are diagonal, I find the following bounds on 

5_i ■ 

the heavy Majorana masses Mj : (a) for a symmetric Dirac neutrino mass matrix 
(no other constraints), a Yb compatible with BBN constraints can be obtained 
for min(M2, M3) > 10 11 GeV (6) if any of the Dirac matrix elements vanishes, 
successful baryogenesis can be effected for a choice of min(M2, M3) as low as a few 
x 10 9 GeV. The latter is compatible with reheat requirements for supersymmetric 
cosmologies with sub-TeV gravitino masses. 



1 Introduction 



The accumulating atmospheric neutrino data from SuperKamiokande jl] has greatly 
increased the likelihood that neutrinos are massive, and that there is mixing among 
the neutrino flavor states. Fits to the zenith angle distribution are consistent with (a) 
maximal v^-v^ mixing and (b) a (mass) 2 difference between the two mass eigenstates 
|Am 2 | ~3x 10~ 3 eV 2 . The solar neutrino data, both from SuperK and other exper- 
iments |], [§], is as yet less definitive in constraining the neutrino masses and mixing: 
there exist the small-angle and large-angle solutions [|] of the MSW effect 0, as well 
as the vacuum oscillation solution. Omitting sterile neutrinos from consideration (per- 
haps with some unwarranted prejudice), one finds that each of these is not yet a clear 
favorite: the day-night asymmetry of solar neutrinos, if persistent at higher statisti- 
cal significance, would disfavor the small-angle MSW and the VO solutions, while the 
reported recoil electron energy spectrum at SuperK requires significant experimental 
or theoretical modification at larger recoil energies in order to be compatible with the 
large-angle solution. It is safe to say that as yet none of the three is ruled out. 

These recent advances have renewed interest in leptogenesis as the precursor to the 
establishment of the cosmological baryon asymmetry Yb = tlbI g*n 1 ~ 0.6 — 1.0 x 10~ 10 
required for a successful description of nucleosynthesis ||. (Here g* is the effective 
number of spin degrees of freedom.) In the simplest leptogenesis scenario @, which 
forms the basis for the discussion in this paper, a B — L asymmetry is established 
through the CP- and L- violating out-of-equilibrium decay of the neutral heavy Ma- 
jorana lepton which partakes in the see-saw mechanism |TIJ for the light neutrino 
masses. (The atmospheric oscillation data indicating very small mass differences can 
be taken as supportive of the see-saw mechanism.) In the next stage of this scenario, 
the B — L asymmetry is reprocessed through the fast (B + L)-violating anomalous 
processes JTT], [H| preceding the electroweak phase transition into the required Yb 
Because the successful completion of this process places non-trivial constraints on the 
both the Dirac and Majorana sectors of the neutrino mass matrix, it has been exten- 



sively discussed in this context in the literature |TJ| . The approach commonly taken is 
to explore the implications for leptogenesis of various models or ansatze for the three 
relevant mass matrices: mp = m^§ ac ,me = m^c rac , and M = M^^ orana . [I is the 
charged lepton.) These constitute important exercises in relating the neutrino data 
and the baryon asymmetry to the underlying Yukawa structures. However, it turns 
out that in certain experimentally allowed regions of the light neutrino mass matrix 
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the analysis becomes much less dependent on the details of the Yukawa matrices, and 
compliance with the BBN-consistent Y# can translate more directly to a constraint 
on the heavy Majorana masses. This is the work of this paper: I will find that for a 
hierarchical light neutrino spectrum, maximal v^-v T mixing, and the small-angle MSW 
solution to the solar neutrino deficit, agreement with the cosmo logical Yb can generally 
obtained for a (lightest) Majorana mass as small as 10 11 GeV. Moreover, when any of 
the matrix elements of mo are zero (in a basis where mi and M are diagonal), this 
bound can be lowered to ~ 10 9 GeV. In the concluding section I discuss the implication 



of these numbers for the gravitino problem [15] of supersymmetric cosmology. 



2 Assumptions 

In what follows, I will work under the following assumptions consistent with present 
neutrino data: 

(a) Light neutrino mass hierarchy |m3| ~5x 10~ 2 eV ^> |m 2 | ^> ni\ 



(b) Maximal v il -v T mixing, consistent with SuperK atmospheric data |L6|, [TJ]] 

(c) Small-angle MSW solution for the solar neutrino deficit 



(d) (Theory) The seesaw mechanism [ID] is operative, with a hierarchical structure 



in the three heavy Majorana masses. 

The hierarchical assumption (a) rules out consideration of a nearly-degenerate scenario 
for neutrino masses in the 2-3 sector. Assumption (c) (the adoption of the small angle 
MSW solution) is not dictated by observation. Compared to the large-angle MSW, it 
provides a marginally better (but not good) fit to the recoil electron energy spectrum 
at SuperK, and a less good fit to the day-night variations ||. More data, and perhaps 
a better understanding of the hep neutrino spectrum, will decide the issue. However, 
it is assumption (c) which allows the study of the constrained system, and that is the 
reason for its adoption |TS| . These assumptions also allow one to ignore renormalization 



group effects in running from the Majorana to low energy scale JT^ . 
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3 Seesaw Relation for Small Angle MSW 



Under Assumption (6) , the electron neutrino plays no role in the atmospheric neutrino 
anomaly, and the light neutrino mixing matrix is given by fl6], p 



U 



( cos 9 sin# \ 

— sin 9/V2 cos9/V2 1/V2 
{ sin6/V2 -cos9/V2 1/V2 j 



where 



U T m u U = diag (e i01 m u e 1 ^ 2 m 2 , m 3 ) . (2) 

The phases (0i,</> 2 ) are Majorana CP-violating phases, and I have omitted the CP- 
violating CKM-type phase in U. The solar mixing angle is given by sin 2 29, and for the 
small-angle MSW solution sin 2 29 ~ 5 x 1(T 3 f|, so that 9 ~ 0.035. 
The seesaw mechanism is expressed by 

m v = rriD M~ x m'^ ) (3) 
= v 2 X M~ l \ T (4) 

where the matrices A and M are defined through the Lagrangian 

C = -LXNH + \NMN , (5) 

so that rriD = Xv, v = 174 GeV. A is a 3 x 3 complex Yukawa matrix. I will work 
in a basis where the charged lepton masses and M are diagonal, so that M~ l = 
diag (1/Mi, I/M2, I/M3). Even if we were to know all six independent elements of 
the symmetric matrix m v as well as the values of Mi, M 2 , M 3 , the seesaw condition (^) 
provides six equations for the nine complex unknown matrix elements of A. In general, 
the leptogenesis scenario requires knowledge of the entire matrix A, so that considerable 
input besides @) is needed in order to determine A. 

The situation improves a great deal for the small-angle MSW solution. I will sim- 
plify matters by taking 9 = in Eq. (HD, and nix — 0, M\ 3> M 2j 3. These are sufficient 
to decouple v e from the z/ M -z/ r seesaw and from the leptogenesis scenario. On input of 
Eqs.(|I]) and (Q), the seesaw equation (f|) (now 2x2) reduces to a set of three equations 
for the four (complex) unknowns a, b, c, d, if we regard m 2 , m^, M 2 , M 3 as 'known': 

(a 2 - c 2 )/M 2 + (d 2 - b 2 )/M 3 = 
a(a- c)/M 2 + d(d-b)/M 3 = m 2 /v 2 

a(a + c)/M 2 + d(d + b)/M 3 = m 3 /v 2 . (6) 
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The Dirac Yukawa matrix A has been parameterized as 

(a d \ 

and the Majorana phase 02 = is incorporated into rri2 : m2 = \m2\e 1 ^. 



4 The Leptogenesis Scenario 

The leptogenesis scenario has been carefully discussed by many authors ||20|| . Briefly, 
the present baryon asymmetry of the universe is calculated in the following manner: 

(a) First, the lepton asymmetry Yl is given in terms of the decay asymmetry e of 
Ni, the lightest of the iV's, parameterized as follows: 

v n L -n L 

Y L = — = k B e , 8 

#*n 7 



where 



1 



M, 



E Irn (AtA)J -± (9) 



167r(AtA)« fr, "-"i v "- vu Mj 



on the assumption of a mass hierarchy Mj M^y^ 0, 21]. The meanings of k 
and B follows: 

(z) Thermal production: If M sma n er = min(M2, Ms) < Trh {Trh= post-inflation 
reheat temperature) and the inverse decay rate is sufficient to establish equi- 
librium, then 



k = suppression factor due to washout by inverse decay 
and 2 — > 2 lepton- violating scattering processes 

B = 1/g* (10) 

where g* is the effective number of massless spin degrees of freedom at the 
time of iV decay (g* = 106.75 in the Standard Model). The factor k is 



determined by numerical integration of the Boltzmann equations || ^] and 
depends most sensitively on the ratio 



K = T t /H(T = M l ) 

\ 8tt ) \ M Pi 
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The suppression factor k reaches its limiting value of 1.0 for K -C 1, and 
drops to ~ 0.01 for K = 20. 

(ii) Non-thermal production via inflaton decay: If M sma ii e r/7W^10 and l^KZilOO, 
then integration of the Boltzmann equations (starting at M/T > 10) reveals 
negligible suppression due to inverse decays, and 



k ~ T RH /M inflaton ~ 10^(T RH /10 iU GeV) 

i? = average number of iV's produced in decay of an inflaton. (12) 
(b) Finally, the baryon asymmetry is established when the B — L asymmetry is 



processed through the fast (B + L)-conserving sphaleron processes |TT|, [12] above 
the electroweak transition temperature, and is given by |23 



B 



\Yl ■ (13) 



5 Results 

It is clear from that a calculation of Yb will involve all the matrix elements of A, 
so that even the truncated 2x2 seesaw equations are not quite sufficient to enable 
a casting of Yg in terms of the masses alone. I will give results for the following 
illustrative constraints on A: 

(1) a, b, c, d, respectively, are set = 0. 

(2) c = d (symmetric case) 

(3) a = b 

Each of these will be worked out for both cases (z) M2 ^> M3 (ii) M3 ^> Mi. In all 
cases, the hierarchy 3> \m,2\ will be respected. 

The results listed in Table I are obtained by inserting these constraints into Eq. (H), 



solving for the matrix elements, and then utilizing Eqs. m) and ( Tip to calculate e and 
the out-of-equilibrium parameter K. The quantities x, y, K , and eo in the Table are 
defined as follows: 

x = \m2\/m^ 

y — -^smaller /-^larger 
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K = m 3 /m , m = (1.7)8iry/g*v 2 /M Pi ~ 1.1 x 1(T 3 eV 

~ 45 

= 3 m 3 M smaller 5 / m 3 \ f M smallcr \ 

60 " 16tt v 2 ~ U.05 eV/UO u GeV/' 1 J 



^smaller — M 3 


^smaller — M 2 




e/eo 


K (^smaller < Trh) 


a = or c = 


b = or d = 


1 


x sin 


4.0 x 10~ 3 


= or d = 


a = or c = 


2x 


| sin 


1.0 x 10~ 2 (x = 0.25) 


c = d or a = b 


c = d or a = 6 


1 


(Axy) 1 ^ 2 sin |0 


4.0 x 10" 3 



Table I: Results for illustrative constraints. For clarity of presen- 
tation, various overall ± signs have been omitted. 



Two comments with respect to the results in Table I are in order: 

• In obtaining the result for e/eo m Line 3, one finds that the algebra simplifies 
greatly if x <C \y. The result given reflects this choice. 

• As noted above, in the case of thermal production, the suppression factor k 
depends on K, and was obtained by integration of the rate equations E 



subject to the initial conditions Y N (M/T ~ 0) = 1/g*, Y L (M/T ~ 0), where 
Y^ = n Ni 1 9* n ^- For non-thermal production, k is given in Eq. (|T2|) above. 

I now proceed to calculate Yb and require 

Y B = F| BN > 0.6 x 10~ 10 . (15) 

From Eqs. (j8|),(|l3]), and (p!4|), one obtains (ignoring signs) 

Y B = 3 x 10- 8 k (Bg*)(e/e ) (^g^) , (16) 

where I have taken m 3 = 5.0 x 10~ 2 eV. From (|TJI) and flT5|), there results a lower 
bound on the lighter of (M 2 , M 3 ): 

k (Bg*) (e/eo) M smaller > 2 x 10 8 GeV . (17) 
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Incorporating the requirements x < \ (mass hierarchy), | sin0| < 1, the three lines 
in Table I can be addressed in turn for each of the scenarios, and a bound obtained 
from Eq. (|TJ: 

Thermal Production: 

Line 1: 

M smaller > 2 x 10 11 GeV . (18) 

Line 2: In this case k depends strongly on x through the dependence of K = 2xK Q , 
and this is reflected in the bound for M sma n cr . 

x = 0.25 : M smaller > 4 x 10 10 GeV 

x = 0.04 : M smaller > 4 x 10 9 GeV (19) 

Line 3: As noted above, the expression for e/eo given in the Table reflects a simplifying 
constraint x \y. The resulting bound is 

M smaller > hy- 1 x 10 10 GeV > 2 x 10 11 GeV (20) 

if we take y < \ in order to maintain the hierarchy in the heavy masses. 

Non-Thermal Production: Here k is given by Eq. (|12|). In terms of a scaling factor 

C = (10 10 GeV/T RH ) (lOOfi)- 1 

one finds 



Line 1: For x < | 



Line 2: 



M smaller > 8 x 10 11 C GeV (21) 



M smaller > 4 x 10 11 C GeV (22) 
Line 3: With the same restrictions as in the previous section, 

M smallcr > 1 x 10 12 C GeV . (23) 
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6 Discussion of Results and Conclusions 



(1) This work has focused on a particular sector of the neutrino mass spectrum (hierar- 
chical) and mixing matrix (maximal Vp-v T mixing, small-angle MSW). With the seesaw 
mechanism, this greatly constrains the Dirac Yukawa matrix A, effectively decoupling 
the electron neutrino and one of the heavy Majoranas. As a consequence, a single 
additional condition on the four complex matrix elements of the effective A allows its 
determination in terms of light and heavy masses, and a single CP- violating Majorana 
phase. The lepton asymmetry resulting from out-of-equilibrium decays of the heavy 
Majorana neutrinos may then be computed in terms of these parameters. Comparing 
the resulting cosmological baryon asymmetry with the value required from BBN then 
turns out to place a lower bound on the lighter Majorana mass. 

(2) In the case of thermal production, for the range of scenarios studied (including 
several not reported in this paper, such as b = c), the lower bounds found for the mass 
of the lightest heavy Majorana are typically of O(10 u ) GeV, well below the inflaton 
mass of ~ 10 13 GeV. Thus, the heavy Majorana may be produced during reheating 
via inflaton decay, without recourse to parametric resonance production |2"5 |. A reheat 
temperature of O(10 u ) GeV requires a large gravitino mass £2.5 TeV [[2f| in order 
that decays of produced gravitinos not destroy the products of nucleosynthesis. If one 
of the entries in the Dirac Yukawa is zero, there are scenarios (Line 2 of Table I), 
in which the smaller Majorana mass may fall below 10 10 GeV, which is a safe reheat 
temperature for low gravitino masses. In the case of non-thermal production, the lower 
bounds are a bit higher, and can exceed the inflaton mass ~ 10 13 GeV if the reheat 
temperature is less than 10 9 GeV (see Eqs. (f21"|)-(|23|)). In that case, production via 
parametric resonance would be necessary. 

(3) Various small parameters, such as the solar mixing angle 9 or the U e s element of 
the mixing matrix, have been set to zero. In principle, small entries for these can 
compete with the mass hierarchy parameters x and y, and cloud the results of this 



work J27J- As a crude measure, one can limit the present discussion to values of 
x,y £ 9 ~ 0.03. However, for the present CHOOZ bound \U e ^\ 2 < 5 x 10~ 2 , a similar 
criterion x, y £ | talmas — 0.2 may be too restrictive. As the data improves, the effects 
of any small non-zero entries can be assessed. 

(4) The discussion presented here is considerably more constrained than previous stud- 
ies jl4| which assume entire textures for both the Dirac and Majorana matrices, and 
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often leave undetermined a good number of parameters. As stated in the introduction, 
such studies are valuable as links to larger theories of flavor symmetries, and are more 
flexible in accommodating a changing scenario for the neutrino parameters. The aim 
here is much more phenomenological, incorporating ab initio certain constraints on 
the light neutrino mass matrix, and leaving to vary only one complex parameter. Of 
course, increased statistics on the day-night effect could begin to seriously disfavor the 
small-angle MSW solution, thus removing the basis for the simplification in this work. 
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